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ABSTRACT 

Using a model of particle acceleration and transport in solar flares, we in- 
vestigate the height distribution of coronal electrons by focusing on the energy- 
dependent pitch-angle scattering. When pitch-angle scattering is not included, 
the peak heights of loop-top electrons are constant, regardless of their energy, 
owing to the continuous acceleration and compression of the electrons via shrink- 
age of magnetic loops. On the other hand, under pitch-angle scattering, the 
electron heights are energy dependent; intermediate energy electrons are at a 
higher altitude, whereas lower and higher energy electrons are at lower altitudes. 
This implies that the intermediate energy electrons are inhibited to follow the 
shrinking field lines to lower altitudes because pitch-angle scattering causes effi- 
cient precipitation of these electrons into the footpoint and their subsequent loss 
from the loop. This result is qualitatively consistent with the position of the 
above-the-loop-top hard X-ray (HXR) source that is located above coronal HXR 
loops emitted by lower energy electrons and microwaves emitted by higher energy 
electrons. Quantitative agreement with observations might be achieved by con- 
sidering primary acceleration before the onset of loop shrinkage and additional 
pitch-angle scattering via wave-particle interactions. 

Subject headings: acceleration of particles — plasmas — Sun: flares — Sun: 
radio radiation — Sun: X-rays, gamma rays 



^Institute for Research on Earth Evolution, Japan Agency for Marine-Earth Science and Technofogy, 
3173-25; Syowa-machi, Kanazawaku, Yokohama 236-0001, Japan 

^ Solar- Terrestrial Environment Laboratory, Nagoya University, Furo-cho, Chikusa-ku, Nagoya 464-8601, 
Japan 



-2- 



Introduction 



One of the most important features of solar flares is the above-the-loop-top 20-50 keV 
hard X-ray (HXR) s ource located higher than coronal < 20 keV HXR loops, discovered by 
Masuda et al.l ( 1l994l ). Although this discovery is regarded as compelling evidence of magnetic 
reconnection, the origin of the above-the-loop-top source itself remains a debatable issue in 
solar flare physics. Following this discovery, the height distribution of coronal HXR sources in 
multiple energy bands has been obser ved and studied using t he Reuven Rama ty Hi gh Energy 



Solar Spectroscopic Imager (RHESST). ISui fc HolmanI (120031 ) . ISui et al.l (120041 ). and lLiu et al. 



( I2OO8I ) discovered two distinct sources near and away from the solar surface. The energy of 
the source near the surface increases toward higher altitudes. This energy-dependent height 
distribution is similar to that of the above-the-loop-top source. The energy of the source away 
from the surface increases toward lower altitudes. Such observations have been regarded as 
evidence of the formation of a current sheet between these two sources. Assuming symmetry 
with respect to the central current sheet, the height distribution of the source away from the 
surface is also similar to that of the above-the-loop-top source. 

Numerous studies on the above-the-loop-top HXR source have been carried out thus 
far. lAschwanden et al.l (jl996al ) analyzed the time-of-fiight differences of HXR light curves, 
and they suggested that the HXR source itself is the electron acceleration site. Plasma 
turbulence, which is potentially generated by a fast flow associated with magnetic reconnec- 
tion, h as been discussed a s a po ssible particle accelerator (IMiller et al.lll996t jPetrosian fc Liu 
20041 ) . iTsuneta fc Naitol ( 1l998l ) suggested that the HXR source is the region where electrons 
are accelerated by fast magnetosonic shocks p roduced by the collisi o n of a fast flow with 
underlying magnetic loops. On the other hand, iMetcalf fc Alexanderl (jl999[ ) analyzed HXR 
spectra at both the corona and footpoint, and they concluded that the HXR sources are a 
consequence of the magnetic trap of parent electrons at the loop top and the subsequent 
pr ecipitation into the footp oint via pitch-angle scattering (termed as trap-plus-precipitation 
bv lMelrose &: Brown (Il976h ). Th e trap-plus-precipitat i on mo del h as been numerical l y stud - 
ied bv Irietcher fc MartensI Jl998[ ). iKarlickv fc Kosugil J2004h . and llVIinoshima et al.l toO^ . 
exhibiting significant confinement of electrons at the loop top (but not necessarily above the 
loop). 

However, previously proposed models have not successfully determined why the HXR 
source at a particular energy band is seen only above the loop and not in the loop, as in the 
case of other energy bands. Recently, iMinoshima et al.l ( 2010l hereafter M2010) developed a 
comprehensive model for particle acceleration and transport in solar flares on the basis of the 
drift-kinetic theory. This model can describe the time evolution of the particle distribution 
function associated with the evolution of flaring electromagnetic fields. As compared to 
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previous models, this model rigorously treats particle acceleration, transport, and dissipation 
processes in a realistic time- varying flare environment. 

To understand the nature of the above-the-loop-top HXR source, we investigate the 
distribution of coronal electrons by using the model of M2010. In particular, we focus on 
the effect of the energy-dependent pitch-angle scattering due to Coulomb collisions on the 
resulting electron height distribution. In Section |2l we briefly describe our model. Simulation 
results are presented in Section [31 We summarize the paper and discuss the nature of the 
above-the-loop-top source in Section HI 



2. Models 



The model solves the drift-kinetic Fokker-Planck equation, 

dt [dt'j di,\dt') dj\dt'J aji I, ""aji, 

where /(r, 7, /i) is the electron distribution function, r is the position in configuration space, 
and (7, fi) are the Lorentz factor and pitch-angle cosine, respectively. The term dr /dt consists 
of the free stream along the magnetic field line and the electric field drift, and the terms d'y/dt 
and dfi/dt consist of the betatron, inertia drift, and magnetic mirror forces (see Equations 
(7)-(9) in M2010). 

In addition, we include the pitch-angle scattering due to Coulomb collisions as a diffu- 
sion term on the right-hand side. W e use the pitch-angle di ffusion coefficient given by the 



Rosenbluth-MacDonald-Judd theory (IRosenbluth et all 119571 ) . 

where K ~ 7.48 x 10"-^^, /3 = (1 —7"^)^'^^, x = {(7 — 1) /T}^^^, and erf(x) is the error 
function. The first term (unity) within the brackets corresponds to the collision with ambient 
cold ions, and the remaining terms correspond to the collisions with electrons themselves. In 
Equation ([2]), it is assumed that the target electrons are almost isotropic in velocity space. 
This is a good approximation because the target electrons are relatively slow and collisional. 
n and T are the electron number density and temperature (average energy) in units of the 
rest mass energy, respectively. They are calculated from the distribution function at each 
time step, 

/I /»oo 
rf/i J d-ff, (3) 
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We use the same model as tha t of electromagn etic fields in M2010. The magnetic field 
model was originally proposed by iLin et al.1 (119951 ): it involves the superposition of two- 
dimensional potential and horizontal fields to impose an X-type neutral line. The simulation 
domain is limited to the area below the neutral line. On the basis of observations (IQiu et al. 
2OO2I ). we change the magnetic field configuration with time to induce electric fields. For 
further details of the model, see M2010. 

In our model, we consider two initial number densities, 10^*^ cm~^ (typical value in active 
regions) and 10^^ cm~^, in order to evaluate the effect of the number density that is used for 
the calculation of the pitch- angle diffusion coefficient. The diffusion term is solved by the 
second-order co nservative form of the interpolated differential operator scheme (IDO-CF; 
Imai et al.ll2008l ) via implicit time integration. Other simulation parameters, the initial and 
boundary conditions, and the numerical method are the same as those in M2010. 



3. Results and Interpretation 

Figured] shows the spatial distribution of the number of 20 keV omnidirectional electrons 
in the cases without (left) and with (right) pitch-angle diffusion (hereafter, the no-diffusion 
and diffusion cases, respectively). The initial number density is 10^° cm~'^. In both the 
cases, the 20 keV electron number is increased around the loop top by betatron acceleration 
and confinement due to shrinkage of the loops, as discussed in M2010. In the diffusion 
case, the electrons are distributed to a greater extent along the field lines toward the lower 
area, as compared to the no-diffusion case. The velocity distribution in the diffusion case 
(not shown) implies that the loss cone is somewhat filled with electrons (depending on 
energy), as compared to that in the no-diffusion case (clear loss-cone distribution. Figure 4 
in M2010). The electron number decreases in the bottom area because electrons in the loss 
cone precipitate into the lower boundary. In the diffusion case, the number further decreases 
because electrons tend to enter the loss cone via pitch-angle scattering. 

One can see that the peak heights of 20 keV electrons at x = (loop top) differ in both 
the cases. The peak height in the diffusion case is higher than that in the no-diffusion case. 
Figure [2] clearly shows the electron height distributions at the loop top. The top and bottom 
plots show the results with initial number densities of 10^" and 10^^ cm~^, respectively. 
In the no-diffusion case (dashed lines), the electron peak heights are constant {z = 0.77), 
regardless of their energy. This can be understood as follows. Any energy electrons are 
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perfectly trapped by the magnetic mirror force without pitch-angle scattering; hence, they 
are continuously accelerated and compressed via loop shrinkage toward a lower altitude. 
Therefore, they peak at the same low altitude. 

On the other hand, in the diffusion case (solid lines), we find that the electron peak 
height depends on the energy of the electrons. Lower energy electrons peak at a lower 
altitude {z < 0.8, close to the no-diffusion case), intermediate energy electrons peak at a 
higher altitude {z > 0.8, far from the no-diffusion case), and higher energy electrons peak at 
a lower altitude {z < 0.8, close to the no-diffusion case). The difference in the peak heights 
is around 1 Mm. The energy of electrons located at a high altitude is higher when the initial 
number density is larger. 

We interpret the simulation r esults in terms of a com bination of loop shrinkage and 



the trap-plus-precipitation model ( iMelrose fc Brown! Il976l ). Unlike the no-diffusion case 



electrons can escape from the magnetic trap via pitch-angle scattering. Some electrons that 
are scattered into a small pitch angle can reach deep into the loop and precipitate into the 
footpoint; finally, they are lost from the loop. The precipitation (loss) rate can be evaluated 
by a combination of the bounce frequency in the loop, Uf, ^ and the pitch-angle diffusion 
frequency, ud ~ -D^^(/i = 0). Here, v is the electron velocity and L is the loop half length. 
When the diffusion frequency is considerably higher/lower than the bounce frequency, the 
precipitation rate Up is approximated by 




for (z/D > (1 - ;Uc) z/fe) , 
for (z/D < (1 - iJ,c) z/fe) , 



where /ic = \/l — Blt/Bfp is the loss-cone angle cosine measured at the loop top, and Blt 
and Bpp are the magnetic field strengths at the loop top and footpoint, resp ectively. Thes e 



two opposite limiting cases are known as the strong and weak diffusion limits f lKennellll969l ). 
Figure E] shows the precipitation rate in the strong (dashed line) and weak (dash-dotted line) 
diffusion limits, and their inverse root mean square (solid line), 

= I (6) 
^(1/z..)^ + {l/u^f 

where n = 10^^ cm~^, L = 15 Mm, and Bfp/Blt = H are the parameters of the loop at 
{x,z) = (0,0.8). 

In the simulation, we consider pitch-angle diffusion via Coulomb collisions. The diffusion 
frequency is considerably higher than the bounce frequency for lower energy electrons; hence, 
their precipitation rate is close to the strong diffusion limit Ug. The strong diffusion limit 
is determined by the electron energy and loop configuration (i.e., the mirror ratio and loop 
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length). In this case, the precipitation rate is lower than the time scale of loop shrinkage 
(dotted line in Figure |3]). The electron loss is low because the precipitation rate is small, 
and the resulting height distribution is similar to that in the no-diffusion case (green lines 
in Figure [2]). Lower energy electrons can follow the shrinking loop to lower altitudes. 

For higher energy electrons, the situation is opposite to that of the lower energy elec- 
trons. The diffusion frequency is considerably lower than the bounce frequency; hence, 
precipitation rate is close to the weak diffusion limit Vyj. The precipitation rate is also lower 
than the time scale of loop shrinkage. In addition, betatron acceleration due to the shrink- 
age increases the number of energetic electrons at the loop top. This contribution is more 
prominent for higher energy electrons because the rate of increase in the electron number 
(9 In // dt) is proportional to the slope of the energy distribution, which is steeper at higher 
energy (for Maxwellian, dlnf/dt oc —dlnf/d'-f = {2K — T)/2KT, where K is the kinetic 
energy). Thus, the electron loss is low, and the resulting height distribution is similar to 
that in the no-diffusion case (pink and red lines in the top plot of Figure [2]). Higher energy 
electrons can also follow the shrinking loop to lower altitudes. 

As shown in Figure [31 the most efficient precipitation takes place in the case of interme- 
diate energy at which the diffusion frequency is comparable to the bounce frequency (a few 
collisions during one loop transit). When the precipitation rate is higher than the time scale 
of loop shrinkage and the rate of increase in the electron number via acceleration, intermedi- 
ate energy electrons are lost from the loop in the course of the shrinkage. This prevents them 
from following the shrinking loop to lower altitudes. As a result, the intermediate energy 
electrons peak at a higher altitude (blue and purple lines in the top and bottom plots of 
Figure ED. 

When the number density is larger than 10^° cm""^, the point of intersection of the dashed 
and dash-dotted lines in FigureOshifts rightward. Thus, the most efficient precipitation takes 
place at higher energy, and this is consistent with the bottom plot of Figure [2], which shows 
that the energy of electrons located at a high altitude is higher (20-50 keV, pink and red 
lines) when the number density is larger (10^^ cm~^). 



4. Summary and Discussion 



We investigated the height distribution of coronal e lectrons in solar flares by using the 
drift-kinetic model developed by lMinoshima et al.l (120101 ). We found that the peak height of 
loop-top electrons depends on their energy under energy- dependent pitch-angle scattering. 
Intermediate energy electrons peak at a higher altitude than higher energy and lower energy 
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electrons because of their efficient loss by precipitation into the footpoint via pitch-angle 
scattering in the course of loop shrinkage. This prevents them from following the shrinking 
field lines to lower altitudes. 

Observations of coronal HXR sources have indicat ed that in the HXR energy range, 
higher energy electrons are located at higher altitudes (IMasuda et al.lll994j : ISui fc Holman 



2003; Sui et al. 2004; Liu et al. 



12008 



Using microwave images taken by the Nobeyama 



Radioheliograph (NoRH). iMasuda et al.l ( l2010l ) statistically investigated the heights of the 
above-the-loop-top HXR and microwave sources. They found that the HXR sources tend 
to be located at higher altitudes than 17 GHz microwaves. Figure S] shows an example 
event, i.e., an X3.1 class flare occurrin g at the we s t limb on 2 002 August 24 (d e tailed reports 
on this flare have been presented by iLi fc GanI (120051 ) and iReznikova et al.l ( 120091 )). It is 
reasonable to assume that electrons producing 17 GHz microw aves are mildly relativistic 
within a typical coronal magnetic field intensity of < 100 Gauss ( lBastiarull999l ). Therefore, 
the observations imply that intermediate energy electrons (several tens to one hundred keV, 
seen as HXR sources) are located at a higher altitude than higher energy electrons (several 
hundred keV, seen as microwave sources) and lower energy electrons (up to several tens of 
keV, seen as coronal HXR loops). Although the evaluation of emissions from the simulated 
electron distribution is necessary, the simulation result including pitch-angle diffusion is 
qualitatively consistent with observations of the energy-dependent height distribution. The 
result without diffusion cannot reproduce the observations because the electrons peak at the 
same altitude, regardless of their energy. 

On the basis of the results, we consider that the above-the-loop-top HXR source is the 
region below which efficient loss takes place for electrons with the highest precipitation rate, 
determined by the balance of the bounce frequency in the loop and the pitch-angle scattering 
frequency. Otherwise, these electrons continue to be trapped in the loop, follow the shrinking 
loop to lower altitudes, and peak at the same altitude observed in other energy bands. 

For numerical simplicity, we do not include the energy change due to Coulomb colli- 
sions in Equation (|T]). The coUisional energy change is important for intermediate energy 
electrons, whereas it is not important for low energy electrons (with energy comparable to 
the temperature of target electrons) and for coUisionless high energy electrons. Interme- 
diate energy electrons may lose energy via collisions. From the viewpoint of particle loss 
and the resulting formation of the height distribution, the following holds in the case of the 
current simulation results: The loss takes place in either velocity or configuration space via 
deceleration or precipitation. 



We have not achieved quantitative agreement between the observations and the model 
thus far. The above-the-loop-top HXR source is typically observed in the 20-50 keV energy 
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band; on the other hand, with an initial number density of 10 cm~ , the current simulation 
has shown that the 10 keV electrons peak at the highest altitude and the 20-50 keV electrons 
peak at lower altitudes (Figure [21 top). This is because betatron acceleration due to loop 
shrinkage significantly increases the number of high energy electrons toward a lower altitude, 
and overcomes the loss via precipitation into the footpoint. For further understanding the 
observations, we consider two possible scenarios. 

In the first scenario, the rate of increase in the high energy electron number via accel- 
eration is reduced if the slope of the energy distribution is fiatter {d\nf/dt cx —dlnf/d'j). 
Thus, electrons should be primarily accelerated before the onset of shrinkage. The initial 
temperature (3 keV) employed in the simulation seems to be insufficient to quantitatively 
explain the observations. If some mechanisms have primarily accelerated the electrons be- 
fore shrinkage, the 20-50 keV electrons can potentially peak at higher altitudes than the 10 
keV electrons because their precipitation rate is not sufficiently small to be approximated by 
the weak diffusion limit, even if it includes only Coulomb collisions (see Figure [3]). Many 
possible scena rios have been p roposed for primary acceleration, such as th e reconnecting 



current sheet (jLitvinenkd Il996l ) and the vicinity of the reconnection region ( iHoshind 12005 
Drake et allbooeh . Evidence of a current sheet has b een observed, and hence, particle ac- 



celeratioii can p otentially occur around this region fISui fc HolmanI l2003l : ISui et all l2004l : 
Liu et al.ll2008l ). iLiu et al.l (|2008[ ) proposed a scenario wherein stronger (weaker) turbulence 
generated at a higher (lower) altitude may result in greater (less) acceleration and energy- 
dependent height distribution of coronal HXR sources. Although it seems reasonable to 
expect stronger turbulence at a higher altitude near the reconnection region, this scenario 
by itself is insufficient for producing microwave-emitting higher energy electrons at a lower 
altitude, where less acceleration is assumed. We believe that a combination of primary ac- 
celeration and subsequent acceleration, transport, and dissipation processes rather than only 
a single acceleration process is crucial to the formation of the above-the-loop-top source. 

In the second scenario, the loss of high energy electrons is enhanced if the precipitation 
rate is increased. As seen in the bottom plot of Figure [21 the 20-50 keV electrons are 
located at higher altitudes than the 10 keV electrons when the initial number density is 
10^^ cm~'^. However, we believe that this value is too high, and hence, unrealistic for a 
typical coronal circumstance. Alternative wave-particle interactions might take place and 
increase the precipitation rate, even if Coulomb collisions are inefficient. 

Let us roughly estimate the pitch-angle scattering rate via wave-particle interactions. 
High-frequency whistler waves are among the most probable agencies of electron scattering. 
The pitch-angle diffusion coefficient due to parallel-propagating whistler waves is approxi- 
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mated by 



TTWce (1 - /i ) k^Wikn) 
D,, ^ (7) 



where Uce is the electron cyclotron frequency, Bq is the backgrou nd magnetic field intensity, 



kn is the resonance wayenum ber, and W{k) is the wave spectrum flKennel fc Petschecklll966 



Kennel Sz Engelmanru Il966l ). The resonance wavenumber satisfies the following dispersion 



relation and cyclotron resonance condition. 



('—] =1 u - fivk 

\U J UJ{U - Uce) 7 



where Up is the plasma frequency and c is the speed of lig ht. Assuming that it has a 



single power-law distr ibution, the wave spectrum is given by ( iHamilton fc Petrosianlll992 



Summers fc Mall2000h 



W{k) = V^It^I V,ot, (9) 



k ■ \ h 



Wtot = / W{k)dk = R-^, 



(10) 



where VFtot is the energy density with a wavenumber greater than the cutoff wavenumber /cmin, 
and R denotes the ratio of the wave energy to the background magnetic energy. Assuming 
that whistler waves are excited via energy transfer from low-frequency Alfven waves and that 
the transfer is much faster than cyclotron damping, we extrapolate the wave spectrum from 
the macroscopic scale to the whistler regime. In the following, we use /cmin = 27r/10^ cm~^ 
and R = lO"! 

In Figure [5|, we compare the pitch-angle diffusion coefficients due to Coulomb colli- 
sions (Equation ([2])) and whistler waves (Equation ([7])) on the basis of different parame- 
ters {uce/ujp, Bo,n, fi,q). Figure [5] (a) shows a reference result with {uce/<^p,Bo,n, fi,q) = 
(0.1,30 Gauss, 9 x 10^ cm""^, —0.7, 1.7). Higher energy ( > 100 keV) electrons are scattered 
by whistler waves rather than by Coulomb collisions; hence, they will be lost to a greater 
extent via precipitation than that expected in the current simulation. This may support the 
second scenario, which suggests that high energy electrons should be lost in the course of loop 
shrinkage. However, the pitch-angle diffusion coefficient due to whistler waves is intricately 
dependent on many plasma parameters. Figure [5] (b)-(d) shows results with higher mag- 
netic field intensity and number density {{BQ,n) = (100 Gauss, 10^^ cm~'^)), steeper wave 
spectrum (g = 2), and higher Uce/^p (0.5) with lower number density {n = 3.5 x 10* cm~'^), 
relative to the reference, respectively. The ratio of the pitch-angle diffusion coefficient due to 
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whistler waves to that due to Coulomb collisions significantly varies according to the plasma 
parameters. It is not easy to evaluate the contribution of wave-particle interactions to the 
total pitch-angle scattering rate in a realistic coronal circumstance. 



Krucker et al.l (120101 ) reported a flare occurring on 2007 December 31, in which a coronal 
microwave source was cospatial with an above-the-loop-top HXR source, as opposed to the 
flare that occurred on 2002 August 24 (Figure H]). On the basis of our model, their result 
may be interpreted as the occurrence of pitch-angle scattering and the resulting efficient 
precipitation into the footpoint at the same altitude (just below the emission peaks) for both 
intermediate and high energy electrons. It is plausible to consider wave-particle interactions 
rather than Coulomb collisions for efficient pitch-angle scattering over a wide range of energy 
(see Figure E]) if sufficient waves are generated. 

The difference in peak heights with different energies is ~ 1 Mm in the simulation, 
whereas the obs erved above-the-loop-top HXR sources are loca ted several Mm higher than 



wnereas tne oDs ervea apove-tne-ioop-top njvti sources are loca tea several Mm nigner tnan 
the flare loops f Masuda et al. 19951 : Aschwanden et al. 1996bl fl). From a combination of 



loop shrinkage and the trap-plus-precipitation model (disregarding the effect of the accel- 
eration), the difference is roughly evaluated as v\Qop [t — {max (z/p)}"^] , where fioop is the 
time-averaged shrinking velocity at the loop top, r is the elapsed time of shrinkage, and 
max{up) is the maximum precipitation rate (~ 0.25 in Figure [3]). If flaring electromagnetic 
flelds are conflgured so as to sustain a fast flow t'loop ~ 1000 km s~^ in r ~ 10 s (time scale 
of nonthermal emissions) or if max (up) is considerably higher, then the difference can be of 
the order of several Mm. If such a condition is not realized, the difference will be < 1 Mm, 
which is difficult to resolve spatially using current instruments. The rare detection of the 
above-the-loop-top source might indicate difficulty in sustaining the fast flow. 

It is noteworthy to suggest that we can observationally determine the maximum precip- 
itation rate if the above-the-loop-top HXR source and flow velocity are observed simultane- 
ously. This is useful for discussing the detailed physical mechanisms of pitch-angle scattering. 
Coulomb collisions, and/or wave-particle interactions. 

On the basis of these discussions, we suggest a formation scenario of the above-the- 
loop-top HXR source, as shown in Figure O For the next step, we consider the acceleration 
near the magnetic reconnection region. Such a study can be combined with the drift-kinetic 
model, which will facilitate a better understanding of the particle acceleration in solar flares. 
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Fig. 1. — Spatial distribution of the number of 20 keV omnidirectional electrons when 
maximum electric fields are generated (5 s after the flare onset, see M2010 for details). The 
left and right images show the results without and with pitch-angle diffusion, respectively. 
The initial number density is 10^° cm~^. The white lines are the magnetic field lines, and 
the dashed lines are the magnetic separatrix. The spatial coordinates are normalized by 
d=15 Mm. 
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Fig. 2. — Electron height distributions at the loop top (x = 0). The solid and dashed lines 
represent the results with and without pitch-angle diffusion. The top and bottom plots show 
the results with initial number densities of 10^° and 10^^ cm~^, respectively. The green, 
blue, purple, pink, and red lines represent electron energies of 1.5, 9.3, 19.5, 48.5, and 79.7 
keV, respectively. The electron number is multiplied by proper factors for illustration. The 
spatial coordinates are normalized by ci = 15 Mm. 
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Fig. 3. — Precipitation rate in the strong (dashed hne) and weak (dash-dotted hne) diffusion 
Umits, and their inverse root mean square (sohd hne, Equation ([H])) as a function of the 
electron energy. A number density of 10^'^ cm~'^ is used. The inverse of the time scale of 
loop shrinkage (5 s)~^ is shown dotted line. 
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Fig. 4.— Spatial distribution of RHESSI 30-50 keV HXR (solid lines) and NoRH 17 GHz 
microwave (dashed lines) overlaid on the Transition Region and Coronal Explorer [TRACER 
195 A negative image of the flare occurring on 2002 August 24. Contour levels are 50%, 
60%, 75%, and 90% of the peak intensity for the HXR, and 20%, 35%, 50%, 75%, and 90% 
for the microwave. 
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Fig. 5. — Comparisons of pitch-angle diffusion coefficients due to Coulomb collisions (Equa- 
tion ([2]), dashed lines) and whistler waves (Equation ([7]), solid lines) as a function of 
the electron energy, (a) Reference result with {uce/<^p, Bq.u, fi,q) = (0.1,30 Gauss, 9 x 
10^ cm~^, —0.7, 1.7). (b)-(d) Results with higher magnetic field intensity and number density 
{{Bo,n) = (100 Gauss, 10^^ cm~^)), steeper wave spectrum (g = 2), and higher Uce/ujp (0.5) 
with lower number density (n = 3.5 x 10® cm~^), relative to the reference, respectively. 
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Fig. 6. — Schematic illustration of the formation scenario of the above-the-loop-top HXR 
source. The solid curved hues denote magnetic field lines, the arrows denote the direction of 
the shrinkage of the loops, the dashed circles denote the emission sources, the cross denotes 
the magnetic reconnection site, and K{Ki < K2 < K3) indicates the electron kinetic energy. 
By efficient precipitation via pitch-angle scattering, electrons with K = K2 do not enter Step 
4 and peak at a higher altitude in the course of loop shrinkage (seen as the above-the-loop-top 
source) than those with K = Ki and K^. 



